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ABSTRACT 

We examine the anisotropies in the power spectrum by the mapping of real to red- 
shift space. Using the Zel'dovich approximation, we obtain an analytic expression for 
the nonlinear redshift space power spectrum in the distant observer limit. For a given 
unbiased galaxy distribution in redshift space, the anisotropies in the power spectrum 
depend on the parameter /(SI) ~ SI 06 , where £1 is the density parameter. We quantify 
these anisotropies by the ratio, R, of the quadrupole to monopole angular moments of 
the power spectrum. In contrast to linear theory, the Zel'dovich approximation predicts 
a decline in R with decreasing scale. This departure from linear theory is due to non- 
linear dynamics and not a result of incoherent random velocities. The rate of decline 
depends strongly on SI and the initial power spectrum. However, we find a universal 
relation between the quantity Rj Ru n (where Ru n the linear theory value of R) and 
the dimensionless variable k/k„i, where k„i is a wavenumber determined by the scale of 
nonlinear structures. The universal relation is in good agreement with a large N-body 
simulation. This universal relation greatly extends the scales over which redshift distor- 
tions can be used as a probe of SI. A preliminary application to the 1.2 Jy IRAS yields 
SI ~ 0.4 if IRAS galaxies are unbiased. 
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INTRODUCTION 



^ Statistical studies of the 

based on galaxy redshift surveys 



large scale structure are mainly 
Redshifts of galaxies dif- 
fer from their distances due to peculiar velocities along the 
line of sight. This difference produces anisotropies in the 
observed galaxy distribution in redshift-space (s-space). The 
anisotropies of clustering in s-space are manifested in the de- 
pendence of the correlation function on the direction of the 
pair separation or, equivalently, in the dependence of the 
power spectrum on the direction of the wavevector. These 
anisotropies offer the promise of measuring the cosmologi- 
cal density parameter, Q, on large scales. This is because 
gravitational instability predicts that the peculiar veloc- 
ity causing the large scale distortions in s-space depend on 
/(Q) = dlnD/dlna « Q° 6 , where D(t) and a(t) are the lin- 
ear growth and scale factors respectively. Therefore, given 
an approximation for gravity, one may obtain an estimate 
of Q by quantifying the anisotropies, say, in the observed 
s-space power spectrum. 

A useful way to characterize the anisotropies in s-space 
is through the angular moments of the s-space power spec- 
trum, P s ; the lowest order moments, the monopole and 
quadrupole, are given by 



+ i 



M 



and 



dfj P s (k, fj) 



+ i 



Q= \ I dn P s (fc, jU )(3 jU 2 - 1) 



(1) 



where /j is the angle between the wavevector k and the line 
of sight. In the linear regime, the ratio R = Q/M is inde- 
pendent of scale and given by 



Ri, 



f/3+f/3 2 
1 + |/3 + ! 



(2) 



where f3 = f(Q)/b, and b is the linear bias factor. Values 
of f3 determined from s-space anisotropies span the range 
~ 0.4- 1.0 (Cole, Fisher, & Weinberg 1995 (CFW); Hamil- 
ton 1995; Fisher, Scharf, & Lahav 1994; Heavens & Taylor 
1995; Fisher et al. 1994; Loveday et al. 1995). One source of 
systematic error in these analyses is in modeling the nonlin- 
ear gravitational evolution (Cole, Fisher, & Weinberg 1994). 
In this Letter, we study the nonlinear evolution of the s- 
space power spectrum. To describe nonlinear dynamics we 
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adopt the Zel'dovich approximation (Zel'dovich 1970). This 
approximation is an excellent tool for describing nonlinear 
dynamics on large scales (e.g., Nusser et al. 1991). The 
Zel'dovich approximation has previously been used by Tay- 
lor (1993) and Schneider & Bartelmann (1995) to study the 
evolution of the real space power spectrum. 

In §2 we derive an analytic expression for the s-space 
power spectrum and find a strong dependence of the ratio R 
on Q and the initial conditions. However, in §3, we show that 
the quantity Rj Ru„ when plotted against the dimensionless 
variable k /k„i , where k„i is a wavenumber determined by 
the scale of nonlinear structures, acquires a universal shape 
independent of Q and the initial power spectrum. This rela- 
tion can be very useful for estimating Q from redshift surveys 
over a wide range of scales. As an example, we present a pre- 
liminary application to the 1.2 Jy IRAS survey. We conclude 
in 84. 



2 THE REDSHIFT SPACE POWER 
SPECTRUM IN THE ZEL'DOVICH 
APPROXIMATION 

Let v, q, x and s be, respectively, the comoving pecu- 
liar velocity, the initial coordinate (i.e, Lagrangian), the 
present real-space (r-space) coordinate (i.e, Eulerian), and 
the s-space coordinate of a particle in units of km s -1 . The 
Zel'dovich approximation states that the displacement vec- 
tor, d, from q to x is related to the peculiar velocity by 
/(Q) - v; accordingly the s-space coordinate, s, is given by 

s = x + (vi)i = q + (I + /(Q)^)d(q,t) , (3) 

where I is the unity matrix and V = 11 T is a matrix which 
projects the velocity vector along the line of sight denoted 
by the unit vector 1. In order to simplify the calculations, we 
work in the "distant-observer" limit where the line of sight is 
well approximated by a fixed direction. Having specified the 
mapping (3), we can readily write down a formal equation 
for the matter density field in s-space as, 



5 S (s) 



d S qS D [s-q-(I + /7>)d]-l 



(4) 



where the s superscript denotes s-space and Sd is the Dirac 
delta function. 

By Fourier transforming (4) we find, 

d 3 s exp[ik • s] S s (s) 

= J d\ exp[ik-q] {exp[ik T (I + /7>)d] - l} (5) 

The second term in brackets reflects the Fourier transform 
of the mean density and contributes only to the k = mode. 
From (5) we see that for k ^ 

^k^k') = J rf 3 qrf 3 q' exp[zk • q+z'k' • q'] 

x <^exp [zk T (I + /7')d + l k' T (I + /7')d']J> , (6) 



where d' = d(q'). The covariance of the modes is related to 
the s-space power spectrum, P s , by 



(S^,) = (2nfS D (k + k')P s (k) , (7) 

where the Dirac delta is a result of translational invariance 
insured by working in the distant observer limit. The rela- 
tions (6) and (7) specify the s-space power spectrum in the 
Zel'dovich approximation. 

The ensemble average in (6) can be written in terms of 
an integral over the joint probability distribution, Pr(d,d'), 
of the displacement vectors 



d 3 dd 3 d'Pr(d,d') 



exp [ik T (I + fV)d + ik' T (I + fV)d' 



(8) 



The great simplification of the Zel'dovich approximation is 
that the displacement is proportional to the initial peculiar 
velocity. If we assume that the initial fluctuations were Gaus- 
sian then, Pr(d,d') is specified by the correlation function 
of the displacement field, C(r) = (dd' T ), where r = q' — q. 
After performing the integral in (8), and comparing with (7) 
we find 



P s (k) 



,3 -ikr 

a re 



x exp[k T (I + /7>)(C(r)-C(0))(I + /7>)k] . (9) 

The displacement field correlation function can be con- 
veniently written in in following form (e.g., Gorski 1988) 

C tJ (r,t) =6 tJ C\(r,t) + r t r 3 [C„ (r, t) - C\(r, t)] ,(10) 

where 



C±(r, t) 
C||(r,t) 



Pit? 
2tt 2 



dk Pu n (k) 



kr 



Jo(kr) 



2ji(kr) 



kr 



where Pu„(k) is the initial r-space power spectrum. 
To first order in C, equation (9) reduces to 



D(t) 2 Pun(k)(l + fii 2 ) 2 



(11) 



(12) 



(13) 



This is the well known result of linear theory (Kaiser 1987); 
the agreement is not surprising since to lowest order the 
Zel'dovich approximation is equivalent to linear theory. 

So far we have assumed that galaxies trace the matter 
distribution. In reality, we should allow for some bias be- 
tween the galaxy and matter distributions. Here we make the 
hypothesis that when the fluctuations are small the galaxy 
and matter density fluctuations are related by a linear bias- 
ing relation. However, with the growth of fluctuations into 
the nonlinear regime at the present time, the bias factor 
becomes a nonlinear function of the density. With this hy- 
pothesis, the parameter /(f2) in (9) is simply replaced by 
Plin = f(Q)/bun, where bu„ is the initial linear bias factor 
scaled to the present using linear theory. 

The integral in (9) can be numerically evaluated to 
yield the nonlinear s-space power spectrum for a given initial 
power spectrum and a value of Q. We introduce a sharp small 
scale cutoff, k c , in the initial power spectrum in an attempt 
to remedy the failure of the Zel'dovich approximation in re- 
gions where orbit mixing has occured. We consider a family 
of CDM-like models specified by the shape parameter T as 
defined in Efstathiou, Bond, & White (1992). We quantify 
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Figure 1. The ratio, R = Q/M, computed from the Zel'dovich 
approximation. The heavy solid and dotted curves correspond to 
an initial power spectrum of (T = 0.5, Q, = 1) with <rg = 1 and 
0.5 respectively. The heavy and light dashed curves correspond 
to (r = 0.2, <T8 = 1.0) with Q, = 1 and 0.3 respectively. The 
truncation of the initial power spectra of these curves is 2tt /k c = 
5 h — 1 Mpc. The light solid curve is the same as the heavy solid 
curve but with truncation at 2ir/k c = 1 h — 1 Mpc. 



the anisotropics in P s at each |k| by the ratio R = Q/M 
where Q and M are computed according to (1). Figure 1 
shows R for a variety of different T, ag, and Q (cf. caption). 
All curves but the light solid one were generated using a 
small scale cutoff of 2n/k c = 5 h -1 Mpc. For comparison, 
the light solid curve was derived with a smaller truncation 
of 2w/k c = 1 h -1 Mpc. We see that, fortunately, R is only 
weakly sensitive to the cutoff. 

On large scales, R asymptotes to Ru„ given in (2), but 
decreases rapidly with increasing wavenumber. This decline 
is because linear theory over predicts the peculiar velocity 
associated with a given density field (Nusser et al. 1991) and 
consequently overestimates the amplitude of anisotropies in 
the power spectrum on scales where nonlinear effects are 
important. Our formalism accounts for multi- valued zones 
in which the flow in r-space is laminar even though particle 
orbits have crossed in s-space. On scales where this effect is 
significant, the quadrupole moment reverses sign and R be- 
comes negative. It is important to distinguish this effect from 
that of incoherent velocities in virialised regions (Fingers of 
God). The scale at which R crosses zero defines a nonlinear 
scale corresponding to a wavenumber k„i. The value of k„i 
is determined by the initial power spectrum as well as Q. 
The Q dependence is due to the fact that for a given r-space 
density field the predicted velocity is proportional to /(O). 
Hence the nonlinear scale in s-space is a function of Q. 



Figure 2. The ratio, R/R\ lrl , from the curves in figure 1 plotted 
against the dimensionless variable r = k/k n \ . The solid dots show 
the empirical fit (14) to the universal relation. The open squares 
with errorbars were obtained from the moments given in CFW 
for their 400 h _1 Mpc, T = 0.25, <r 8 = 0.8, and V. = 0.3 N-body 
simulation (cf., CFW figure 2). 



3 UNIVERSAL SCALING RELATION 

In the Zel'dovich approximation the peculiar velocity is 
specified by the r-space density field independent of the ini- 
tial conditions (Nusser et al. 1991). Since the anisotropies 
in s-space are a reflection of this density /velocity relation, 
we expect the ratio R to be a weak function of the initial 
power spectrum when expressed in terms of the dimension- 
less variable r = k/k„i. In order to verify this ansatz, we 
plot, in figure 2, the ratio R/Ri, n against the variable r 
for each of the curves in figure 1. In addition, we show re- 
sults from the large N-body simulation of CFW. This figure 
shows that, indeed, the ratio R/Ri, n acquires a universal 
shape over a large range of scales (over two orders of mag- 
nitude). The agreement of the Zel'dovich predictions with 
the simulation is remarkable. The small differences between 
the different curves may be attributed to two effects. First, 
the nonlinear density /velocity relation in Fourier space in- 
volves coupling between different modes; this coupling may 
depend on the initial power spectrum. Second, the failure 
of Zel'dovich approximation in regions where orbit mixing 
has occured is especially pronounced for models with large 
power on small scales; in reality, however, caustics are pre- 
served and therefore one might expect the scaling relation 
to hold more tightly. 

A simple empirical fit to the scaling relation, shown by 
the solid dots in figure 2, is 



R 

Rim 



(14) 



1 + fr 3 / 2 

This relation permits an estimation of the parameter Q or 
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0.5 

T 



km s -1 . Yet, rich clusters are rare in the universe and can, in 
principle, be excised from the analysis. Field galaxies have 
a random velocity component in addition to their coher- 
ent motions. However, the amplitude of this component is 
small. The one-dimensional velocity dispersion of all galax- 
ies, including those in clusters, is estimated to be PS 140-200 
km s" 1 (Miller, Davis & White 1995). 

Measurements of the redshift distortions in the planned 
Sloan Digital Sky Survey (Gunn & Knapp 1993) and Anglo- 
Australian 2dF galaxy survey should provide accurate deter- 
minations of R over a wide range of scales. Our formalism 
is especially suitable for the estimation of Q from these sur- 
veys. 
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Figure 3. The ratio, R, from the 1.2 Jy IRAS survey plotted 
against the variable r. The points were obtained from figure 7 
in CFW. The curves correspond to the universal scaling model 
(equation 14) with values of pi tn = 0.3, 0.6, and 1.0 as denoted 
on the figure. 

more precisely f3u„ on scales where nonlinear effects are im- 
portant. Here we use the relation (14) to estimate f3u„ from 
the 1.2 Jy IRAS survey (Fisher et al. 1995). The ratio R 
for this survey versus scale was taken from CFW (figure 7). 
The wavenumber at which R crosses zero was estimated to 
be k n i = 0.22 hMpc -1 . The three solid curves in the figure 
show the value of R determined from the scaling relation 
(14) for different f3u„. Though the errors in R are large, the 
data favor f3u„ ~ 0.6. 



4 DISCUSSION 

We have presented a calculation of the anisotropies in the 
s-space power spectrum based on the Zel'dovich approxima- 
tion. The Zel'dovich approximation shows clearly that the 
decline of the quadrupole to monopole ratio on large scales 
is due to nonlinear, yet coherent, peculiar velocities. Despite 
a strong dependence on Q and the initial power spectra, the 
s-space anisotropies obey a universal scaling relation when 
expressed as a function of the variable k /k„i . This scaling 
relation extends analyses of redshift distortions into the non- 
linear regime with the only free parameter being f3u„. This 
is a marked improvement over phenomenological models for 
the distortion which require an unrealistically high random 
velocity component to explain the decline of R with scale 
(e.g., Peacock & Dodds 1994; CFW). 

We have not incorporated the effects of random veloci- 
ties. The main manifestation of random velocities is the so- 
called Finger of God effect arising from the redshift stretch- 
ing of virialised regions of rich clusters. The velocity dis- 
persion in these regions can be on the order of ~ 1000 
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